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Abstract 

We present a new leading and next-to-leading QCD analysis of the world data 
on inclusive polarized deep inelastic scattering. A new set of polarized parton 
densities is extracted from the data and the sensitivity of the results to the newly 
incorporated SLAC/E155 proton data is discussed. 



PACS numbers:13.60.Hb; 13.88-he; 14.20.Dh; 12.38.-t 



1 Introduction 



Deep inelastic scattering (DIS) of leptons on nucleons has remained the prime source 
of our understanding of the internal partonic structure of the nucleon and one of the 
key areas for the testing of perturbative QCD. Decades of experiments on unpolarized 
targets have led to a rather precise determination of the unpolarized parton densities. 
Spurred on by the famous EMC experiment ^ at CERN in 1988, there has been a huge 
growth of interest in polarized DIS experiments which yield more refined information 
about the partonic structure of the nucleon, i.e., how the nucleon spin is divided up 
among its constituents, quarks and gluons. Many experiments have been carried out 
at SLAG, CERN and DESY. 

In this paper we present an updated version of our NLO polarized parton densities 
in both the MS and the JET (or so-called chirally invariant) factorization schemes 
as well as the LO ones determined from the world data ^, Q on inclusive polarized 
DIS. Comparing to our previous analysis P: 

i) For the axial charges 03 and Og their updated values are used |p, 0]: 

ag = = F + D = 1.2670 ± 0.0035, ag = 3F - D = 0.585 ± 0.025 . (1) 

ii) In our ansatz for the input polarized parton densities 

Af,{x,Ql)=p,ix) /f^^^(x,gg) (2) 

we now utilize the MRST'99 set of unpolarized parton densities fi{x, Ql) instead of 
the MRST'98 one. 

iii) The recent SLAC/E155 proton data are incorporated in the analysis, 
iv) The positivity constraints on the polarized parton densities are discussed. 



2 Method of Analysis 

The nucleon spin-dependent structure function of interest, g^{x,Q^), is a linear com- 
bination of the asymmetries and (or the related virtual photon-nucleon asym- 
metries ^^2) measured with the target polarized longitudinally or perpendicular to 
the lepton beam, respectively. Neglecting as usual the subdominant contributions, 
A^{x,Q^) can be expressed via the polarized structure function g^{x,Q^) as 

<(^.«^)-(i + 7^)||r|; (3) 
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or 

Q') = jw^^M^ + R^'i^^ Q')] (4) 

using the relation between the unpolarized structure function Fi{x, Q"^) and the usually 
extracted from unpolarized DIS experiments F2{x,Q'^) and R{x,Q^) 

2xF^ = F^{1 + 7')/(l + R"^) {N = p, n, d). (5) 

In Eq. the kinematic factor 7^ is given by 

7^^^^. (6) 

In Mn is the nucleon mass. It should be noted that in the SLAG and HERMES 
kinematic region 7 cannot be neglected. 

In the NLO QCD approximation the quark-parton decomposition of the proton 
structure function g\{x, Q^) has the following form (a similar formula holds for (7"): 

Q') = \i: el[{^<l + ® (1 + ^^5C,) + AG ® (7) 

Z g Ztt Ztt l\f 

where Ag(x, Q^), Ag(x, Q^) and AG{x,Q^) are quark, anti-quark and gluon polar- 
ized densities in the proton which evolve in according to the spin-dependent NLO 
DGLAP equations. SCq^c are the NLO terms in the spin-dependent Wilson coefficient 
functions and the symbol denotes the usual convolution in Bjorken x space. Nf is 
the number of flavours. 

It is well known that at NLO and beyond, the parton densities as well as the 
Wilson coefficient functions become dependent on the renormalization (or factorization) 
scheme employed.^] Both the NLO polarized coefficient functions and the NLO 
polarized splitting functions (anomalous dimensions) |]10| needed for the calculation of 



gi{x,Q'^) in the MS scheme are well known at present. The corresponding expressions 
for these quantities in the JET scheme can be found in ||11|| . 



All details of our approach to the fit of the data are given in ||T2[. Here we would 

like to emphasize only that according to this approach flrst used in fl^, the NLO 

QCD predictions have been confronted to the data on the spin asymmetry A^{x, Q"^), 

rather than on the gi{x, Q"^). The choice of appears to minimize the higher twist 

(HT) contributions to which are expected to partly cancel with those of m 

the ratio (^, allowing use of data at lower (in polarized DIS most of the small 

^Of course, physical quantities such as the virtual photon- nucleon asymmetry Ai{x,Q^) and the 
polarized structure function gi{x,Q^) are independent of choice of the factorization convention. 
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X experimental data points are at low Q"^). Indeed, we have found |0 that if for gi 
and Fi leading-twist (LT) QCD expressions are used in (j^), the HT corrections h{x) to 
= Ai{x,Q'^)ir[ + h{x)/Q'^, extracted from the data, are negligible and consis- 
tent with zero within the errors (see Fig. 1). (Note that the polarized parton densities 
in QCD are related to the leading- twist expression oi gi-) On other hand, it was shown 
||T5| that if F2 and R in Eq. (H) are taken from experiment (as has been done in some 
of the analyses) the higher twist corrections to Ai are sizeable and important. So, in 
order to extract the polarized parton densities from gi data the HT contribution to gi 
(unknown at present) has to be included into data fit. Note that a QCD fit to the gi 
data keeping in gi{x, Q^)qcd only the leading-twist expression leads to some "effective" 
parton densities which involve in themselves the HT effects and therefore, are not quite 
correct. These results suggest that in order to determine polarized parton densities less 
sensitive to higher twist effects, it is preferable at present to analyze Ai data directly 
using for gi and Fi their leading twist expressions. Bearing in mind this discussion one 
must be careful when comparing the polarized parton densities determined from the 



inclusive DIS data with those extracted from semi-inclusive DIS data 161. The results 



will depend, especially in a LO QCD analysis, upon whether the higher twist terms 
are or are not taken into account. 



Following the procedure of our previous analyses we have extracted the NLO (as 
well as LO) polarized parton densities from the fit to the world data on A^{x,Q'^) 
using for the flavour non-singlet combinations of their first moments 

03 = (Am + Am) (g') - (Ad + M) (g") , (8) 

as = (Am + An)(g') + (Ad + M){Q^) - 2(As + As)(g2) , (9) 

the sum rule values ([l|). The sensitivity of the polarized parton densities to the devia- 
tion of flg from its SU(3) flavour symmetric value (0.58) has been studied and the results 
are given in Here we will present only the polarized parton densities corresponding 
to the SU(3) symmetric value of ag in (|l]). 

What we can deduce from inclusive DIS in the absence of charged current neutrino 
data is the sum of the polarized quark and anti-quark densities 

(Am + AM)(a;,g2), (Ad + A(J)(a;, g^), (As + As)(a;, g^) (10) 

and the polarized gluon density AG'(x, g^). The non-strange polarized sea-quark den- 
sities A'u(x,g^) and A(i(x,g^), as well as the valence quark densities AM^(a;, g^) and 
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Au^ = Au- Au, Ad^ = Ad- Ad (11) 

cannot be determined without additional assumptions about the flavour decomposition 
of the sea. Nonetheless (because of the universality of the parton densities) they are of 
interest for predicting the behaviour of other processes, like polarized pp reactions, etc. 
That is why, we extract from the data not only the quark densities (0) and AG'(x, Q^), 
but also the valence parts Au„(a;,(5^), A(i„(x,Q^) and anti-quark densities using the 
assumption on the flavour symmetric sea 

AUsea = Au = Adsea = Ad = As = As. (12) 

For the input LO and NLO polarized parton densities at Qo = 1 GeV^ we have 
adopted a very simple parametrization 

xAdv{x,Ql) = r]dAdx''''xdv{x,Ql), 



xA,s(x, Qq) — ffgAgX "xs^XjQq), 

xAG{x,Ql) = r]gAgx''^xG{x,Ql), (13) 
where on RHS of (|13D we have used the MRST98 (central gluon) pl and MRST99 



(central gluon) P] parametrizations for the LO and NLO (MS) unpolarized densities, 
respectively. The normalization factors ^/ in (|13D are fixed such that rjf are the first 
moments of the polarized densities. The first moments of the valence quark densities 
T]u and rjd are constrained by the sum rules (|I]) . The rest of the free parameters in (^3|) , 

{a„, ad, r]s, a^, T]g, a^}, (14) 

have been determined from the best fit to the data. 

The guiding arguments to choose for the input polarized parton densities the ansatz 
([T3| ) are simplicity (not too many free parameters) and the expectation that polarized 
and unpolarized densities have similar behaviour at large x. Also, such an ansatz allows 
easy control of the positivity condition, which in LO QCD implies: 

|A/,(x,Q^)| < f^{x,Ql), \AMx,Ql)\ < Mx,Ql). (15) 



The constraints (p^ ) are consequence of a probabilistic interpretation of the par- 
ton densities in the naive parton model, which is still correct in LO QCD. Beyond 
LO the parton densities are not physical quantities and the positivity constraints on 
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the polarized parton densities are more complicated. They follow from the positivity 
condition for the polarized lepton-hadron cross-sections Ao"j with the unpolarized ones 
(lAcTil < cTj) and include also the Wilson coefficient functions. It was shown how- 
ever, that for all practical purposes it is enough at the present stage to consider LO 
positivity bounds, since NLO corrections are only relevant at the level of accuracy of 
a few percent. 

3 Results 

In this section we present the numerical results of our fits to the world data on 
y4i(x, Q^). The data used (185 experimental points) cover the following kinematic 
region: 

0.005 <x< 0.75, 1 < < 58 GeV^ . (16) 

The total (statistical and systematic) errors are taken into account. The systematic 
errors are added quadratically. 

We prefer to discuss the NLO results of analysis in the JET scheme. To compare our 
NLO polarized parton densities with those extracted by other groups, we present them 
also in the usually used MS scheme using the renormalization group transformation 
rules for the parton densities. 

It is useful to recall the transformation rules relating the first moments of the 
singlet quark density, AS((5^), and the strange sea, (As + As){Q'^), in the JET and 
MS schemes: 

ASjET = A^MsiQ') + Nj^^AGiQ'), (17) 

{As + As)jET = {As + As)j^{Q') + ^f^AG{Q'), (18) 

where AG{Q'^) is the first moment of the polarized gluon density AG{x, Q^) (note that 
AG is the same in the factorization schemes under consideration). 

A remarkable property of the JET (and so-called Adler-Bardeen (AB) [|T9|) schemes 
is that the singlet AT,{Q'^), as well as the strange sea polarization (As + As){Q'^), are 

independent qusmtities. Then, in these schemes it is meaningful to directly interpret 
AS as the contribution of the quark spins to the nucleon spin and to compare its value 
obtained from DIS region with the predictions of the different (constituent, chiral, etc.) 
quark models at low Q^. 

It is important to mention that the difference between the values of the strange 
sea polarization, obtained in the MS and JET schemes could be large if AG in ([l8|) is 
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positive and large. To illustrate how large it can be, we present the values of (As + As) 
at = 1 GeV^ obtained in our analysis of the world DIS data in the MS and JET 
schemes (AG = 0.68): 

(As + As)ms = -0.13 ±0.04 

(As + As)jET = -0.07 ±0.02. (19) 

Note that if AG is larger than 0.68, (As ± As) jet could vanish in agreement with what 
is intuitively expected in quark models at low-Q^ region (Q^ ~ 0). 

The numerical results of our fits to the data are summarized in Table 1. The best 
LO and NLO(JET) fits correspond to per degree of freedom of Xdflo — 0.921 and 
to Xdfnlo — 0.871. In LO QCD AG{x,Q'^) does not contribute directly in gi and the 
gluons cannot be determined from DIS data alone. For that reason the LO fit to the 
data was performed using for the input polarized gluon density AG(a;, Ql) that one 
extracted by the NLO fit to the data: 

AG(x, QDlo = AG(x, g^)NLO(jET) . (20) 

We consider that such a procedure leads to non-realistic errors of the rest of parameters 
and therefore, we present only their central values in Table 1. It is important to note 
that in the polarized case the LO approximation has some peculiarities compared to 
the unpolarized one. As has been already mentioned above, as a consequence of the 
axial anomaly, the difference between NLO anti-quark polarizations Ag^ in different 
factorization schemes could be quite large, in order of magnitude of Ag^ themselves. 
In this case the leading order will be a bad approximation, at least for the polarized 
sea-quark densities extracted. Also, bearing in mind that in polarized DIS most of the 
data points are at low Q^, lower than the usual cuts in the analyses of unpolarized 
data {Q^ > 4 — 5 GeV'^), the NLO corrections to all polarized parton densities are 
large in this region and it is better to take them into account. Nevertheless, the LO 
polarized parton densities may be useful for many practical purposes: For preliminary 
estimations of the cross sections in future polarized experiments, etc. For that reason 
we present them in this paper, at = 1 GeV^ (see Fig. 2 and Appendix) and for any 

in the region: 1 < < 6.10^ (the FORTRAN code is available^) . Further, we will 
discuss mainly the NLO QCD results. 

As in our previous analysis a very good description of the world data on 
and Qi is achieved. The new NLO theoretical curves for Ai corresponding to the 

■fhttp://durpdg.dur.ac.uk/HEPDATA/PDF 
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best fit practically coincide with the old ones (see Fig. 3). The agreement with the 
SLAC/E155p data involved in this analysis is also very good. This is illustrated in Fig. 
4. The NLO QCD theoretical predictions for Ai corresponding to our previous analysis 
of the data are also shown. They are in a good agreement with the SLAC/E155 data 
not available at the time our previous analysis was performed. In Fig. 5 we compare the 
new NLO gi curves with the SMC and SLAC/E143 proton data. The extrapolations 
for Qi in the yet unmeasured small x region at different are also shown. As seen from 
Fig. 5, the proton spin dependent structure function changes sign at x smaller than 
10^^ and becomes negative if the gluon polarization is positive. One of the challenges to 
future polarized DIS experiments is to confirm or reject this behaviour, quite different 
from the usual Regge type behaviour. 

The extracted NLO (JET) polarized parton densities at Q"^ = 1 GeV^ are shown in 
Figs. 6(a) and 6(b). (The explicit expressions are given in the Appendix.) The new 
parton densities are found to be within the error bands of the old ones. The positivity 
constraints have not been imposed during the fit. Except for the strange sea density 
As{x), the polarized parton densities determined from the data are compatible with 
the LO positivity bounds ([ISD imposed by the MRST99 unpolarized parton densities. 
However, if one uses the more accurate LO positivity bounds on As{x) obtained by 
using the unpolarized strange sea density s{x)bpz (Barone et al. 0]), As(x) also lies 
in the allowed region. It is important to mention that s{x)bpz is determined with a 
higher accuracy compared to other global fits. 

In Fig. 7 we illustrate the difference between the singlet quark density AS(x,(5^) 
and the strange sea density (As+As)(x, Q^) determined in the JET and MS schemes at 
Q"^ = 1 GeV^. They differ essentially in the small x region and this difference increases 
with AG increasing. In Figs. 8(a) and 8(b) we compare our NLO (MS) polarized parton 
densities with those obtained by AAC |0| and GRSV using almost the same set 
of data. This comparison is a good illustration of the present situation in polarized 
DIS. While the quality and the kinematic range of the data are sufficient to determine 
the valence polarized densities Am^,(x, Q^) and AMt,(x, Q^) with a good accuracy (if 
an SU(3) symmetry of the flavour decomposition of the sea is assumed), the polarized 
strange quark density As(x, Q^) as well as the polarized gluon density AG(x, Q^) are 
still weakly constrained, especially AG. 

Finally, let us turn to the quark and gluon polarizations (the first moments of the 
polarized parton densities). The results of the new as well as of the old analysis are 
presented in Table 2 (Q^ = 1 GeV^). The corresponding values in the MS scheme are 
also given. The changes of the central values of the parton polarizations (JET scheme) 
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are negligible and within the errors of the quantities. Note that for the central value of 
the axial charge ao{Q^) (equal to AE(Q^)^ in MS scheme) we obtain now somewhat 
smaller value: ao(l GeV^) = 0.21±0.10 than the old one: ao = 0.26±0.10. The values 
of the LO parton polarizations are between those of the JET and MS schemes. 

For the gluon polarization AG corresponding to the best NLO(JET) fit we have 
found AG = 0.68 ± 0.32 at = 1 GeV^. However, if one takes into account the sen- 
sitivity of AG to variation of the non-singlet axial charge Og from its SU(3) symmetric 
value of 3F-D, the positive values of AG could lie in the wider range [0, 1.5] |T^. A 
negative AG is still not excluded from the present DIS inclusive data. 

4 Conclusion 

We have re-analyzed the world data on inclusive polarized deep inelastic lepton-nucleon 
scattering in leading and next-to-leading order of QCD adding to the old set of data 
the SLAC/E155 proton data. Compared to our previous analysis: i) the updated 
values for the non-singlet axial charges qa and ag and ii) in the ansatz for the input 
polarized parton densities the MRST99 set of unpolarized parton densities instead of 
the MRST98 one have been used. It was demonstrated that the polarized DIS data 
are in an excellent agreement with the pQCD predictions for Ai{x, Q^) and gi{x, Q^) 
and that the new theoretical curves practically coincide with the old ones. A new set 
of NLO polarized parton densities in the JET and MS factorization schemes as well 
as LO polarized patron densities have been extracted from the data. We note that 
one must be careful when using LO polarized densities. The LO QCD approximation 
could be a bad one, at least for the polarized sea-quark densities, if it will turn out 
that the gluon polarization is positive and large (the first direct measurement of AG 
p2| as well as the QCD analyses support such a possibility). We have found that the 
new NLO (JET) polarized parton densities do not change essentially and lie within the 
error bands of the old parton densities. 

What follows from our analysis is that the limited kinematic range and the pre- 
cision of the present generation of inclusive DIS experiments are enough to deter- 
mine with a good accuracy only the polarized parton densities {Au + Au){x, Q"^) and 
{Ad + Ad){x, Q^). The polarized strange sea density (As + As){x, Q^) as well as the 
polarized gluon density AG(x, Q"^) are still weakly constrained, especially AG. The 
non-strange polarized sea-quark densities Au and Ad cannot be determined, in prin- 
ciple, from the inclusive DIS experiments alone without additional assumptions. The 
further study of flavour decomposition of the sea as well as a more accurate determi- 
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nation of the gluon polarization are important next steps in our understanding of the 
partonic structure of the nucleon and this will be done in the forthcoming and future 
polarized lepton-hadron and hadron-hadron experiments. 



A FORTRAN package containing our NLO polarized parton densities in both the 
JET and MS factorization schemes as well as LO parton densities can be found at 
[http: / / durpdg.dur.ac.uk/HEPDAlA / PDF] or obtained by electronic mail: 



sidorov@thsunl.jinr.ru or stamenov@inrne.bas.bg. 
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Appendix 

For practical purposes we present here explicitly our LO and NLO (JET) polarized 
parton densities at = 1 GeV^. The polarized valence quark densities correspond to 
SU(3) flavour symmetric sea. 

LSS - LO: 

xAu^{x) = 0.3112 (1 - x)^-^^^ ( 1 - 0.4085 x^^^ + 17.60 x ) , 

xAd,,{x) = -0.01563 x"-^^^" (1 - xf^^^ ( 1 + 37.25 x^^^ + 31.14 x ) , 

xAs{x) = -0.08548 (1 - xf^^'"^ ( 1 - 0.9052 x^^^ + 11.53 x ) , 

xAG{x) = 19.14 x^-^^^ (1-x)^-^^^ ( 1-3.147x^/2 + 3.148 0; ) . (A.l) 



LSS - NLO(JET): 

xAu^{x) = 0.5052 (1 - x)^-^^^ ( 1 + 2.179 x^/^ + 14.57 x ) , 

xAd^ix) = -0.01852 x°-27°^ (1 - x)^-^^^ ( 1 + 35.47 x^/^ + 28.97 x ) , 

xAs(x) = -0.1525 x^-^^2 (1 - x)^-^^^ ( 1 + 3.656 x^/^ + 19.50 x ) , 

xAG'(x) = 19.14x1-"^ (1-x)^-®^^ ( 1-3.147 xi/2 + 3.148x ) . (A.2) 
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Table 1. Results of the LO and NLO fits to the world data {Ql = 1 GeV^). 
The errors shown are total (statistical and systematic) . The parameters marked 
by (*) are fixed. 



Fit 


LO 


NLO (JET) 


DF 


185 - 4 


185 - 6 




166.7 


155.9 




0.921 


0.871 




0.926* 


0.926* 




0.121 


0.253 ± 0.027 


Vd 


- 0.341* 


-0.341* 


ad 


0.102 


0.000 ± 0.054 


Vs 


- 0.055 


- 0.036 ± 0.007 




0.754 


1.613 ± 0.429 


Vg 


0.681* 


0.681 ± 0.141 


ag 


0.149* 


0.149 ± 0.741 



Table 2. First moments (polarizations) of polarized parton densities at = 1 GeV'^. 



Fit 


Aii + Am 


M + M 


As-F As 


AG 


AE 


Oo 


okl/JET 


0.8(3±0.03 


-0.10±0.05 


-0.06±0.02 


0.57±0.31 


0.10±0.06 


0.26±0.10 


new/JET 


0.85±0.03 


-0.41±0.05 


-0.07±0.02 


0.68±0.32 


0.37±0.07 


0.21±0.10 


new/ MS 


0.80 


-0.47 


-0.13 


0.68 


0.21 


0.21 


LO 


0.82 


-0.45 


-0.11 


0.68 


0.25 


0.25 
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Figure Captions 

Fig. 1. Higher twist contribution [x) to tlie spin asymmetry Q^) extracted 

from tlie data. 

Fig. 2. Leading order polarized parton densities at = 1 GeV^. 

Fig. 3. Comparison of our NLO results in JET scheme for A^{x^Q'^) with the 
experimental data at the measured x and values. Errors bars represent the total 
(statistical and systematic) error. Our old NLO results are shown for comparison. 

Fig. 4. Comparison of our NLO(JET) result for A\ with SLAC/E155p experimental 
data. Error bars represent the total errors. The predictions for A\ (dot curves) from 
our old analysis |^ are also shown. 

Fig. 5. Comparison of our NLO(JET) resuhs for gl with SMC and SLAC/E143 pro- 
ton data. The extrapolations at small x are also shown. 

Fig. 6(a). NLO(JET) polarized parton densities x(Am + Au) and x{Ad + Ad) at 
Q"^ = 1 GeV^. The old parton densities together with their error bands are presented 
for comparison. 

Fig. 6(b). NLO(JET) polarized strange sea x{As + As) and gluon polarized parton 
densities at = 1 GeV^. The old parton densities together with their error bands are 
presented for comparison. 

Fig. 7. Comparison between NLO polarized singlet and strange sea parton densities 
at g2 = 1 GeV^ in the JET and MS schemes. 

Fig. 8(a). Comparison between our NLO(MS) polarized valence quark densities at 
= 1 GeV2 with those obtained by AAC (NLO-2 set) [0 and GRSV ('standard' 



scenario) [|T^ 



Fig. 8(b). Comparison between our NLO(MS) polarized strange quark and gluon 
densities at = i QeV^ with those obtained by AAC (NLO-2 set) M and GRSV 



[^'standard' scenario) |T^ 
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